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FOREWORD 


This  report  was  prepared  under  Task  ZR0000101,  IR-59  and  OMN  funds. 

It  is  concerned  with  the  relationship  between  infinite  diameter  detonation 
velocity  and  loading  density  of  high  explosives  and  mixtures  of  high 
explosives.  Methods  for  predicting  the  velocity  of  the  latter  are  also 
considered.  The  results  should  be  of  interest  to  the  areas  of  detonation 
chemistry  and  physics;  they  are  also  relevant  to  the  fields  of  safety  and 
sensitivity. 
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THE  DETONATION  VELOCITY-LOADING  DENSITY  RELATION 
FOR  SELECTED  EXPLOSIVES  AND  MIXTURES  OF  EXPLOSIVES 

INTRODUCTION 


For  many  years.  It  has  been  accepted  that  detonation  velocity  (D)  is  a 
linear  function  of  the  loading  density  (p0)  In  the  range  p0  2  1.  Indeed,  In 
1945,  Hurwitz1  published  linear  relqtlons  derived  from  the  extensive  ERL  data 
for  numerous  explosives.  His  straight  lines  were  drawn  through  data  of  charges 
prepared  from  the  finest  particle  size  HE  and  with  the  largest  diameter  at 
p0  *  0.8  g/cm3.  Hurwltz's  values  of  slopes  and  intercepts  have  withstood  the 
test  of  time  very  well;  they  are  quoted  In  at  least  two  texts**3  and  In  many 
reports.  However  in  1961,  It  was  reported  that  the  D  vs  p0  curve  of  TNT4 
showed  a  sharp  change  in  slope  at  po  *  1.5324  g/cm3.  About  nine  years  later, 
a  similar  change  In  slope  was  reported  In  the  curve  for  PETN5  at  p0  •  1.65 
g/cm3. 


There  Is  a  possible  explanation  for  this  deviant  behavior,  and  one 
objective  of  this  work  Is  to  point  It  out.  A  second  Is  to  find  an  acceptable 
way  to  predict  the  detonation  velocity  of  a  porous  mixture  from  a  knowledge 
of  0(po)  of  the  components.  Finally,  the  predictions  will  be  compared  with 
experimental  values  and  with  approximations  obtained  by  several  methods. 

In  efforts  to  obtain  high  density  pressed  charges,  some  Investigators 
have  used  heated  molds  or  solvents  or  both  In  the  pressing  procedure.  In 
general,  there  is  no  record  of  how  this  would  make  the  charge  behavior  differ 
from  that  of  the  more  widely  used  cold-pressed  charges.  In  particular.  It 
nv’.y  Increase  the  critical  diameter  dc  and  hence  the  diameter  effect  on  D  of 
any  size  charge.  These  effects  would  be  expected  for  HE  such  as  TNT  and  PETN 


^Hurwitz,  M.  D.,  "The  Rate  of  Detonation  of  Various  Explosive  Compounds 
and  Mixtures,"  OSRD  5611  (1945). 

2Cook,  M.  A.,  The  Science  of  High  Explosives  (New  York:  Reinhold,  1958). 
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Johansson,  C.  H. ,  and  Persson,  P.  A.,  Detonics  of  High  Explosives 
(New  York:  Academic  Press,  1970). 

*Ur1zar,  M.  J. ,  James,  Jr.,  E. ,  and  Smith,  l.  C.,  "Detonation  Velocity  of 
Pressed  TNT,"  Phvs.  Fluids,  Voi.  4,  No.  2,  1961,  pp. 262-274. 

^Hornlg,  H.  C. ,  Lee,  E.  L.,  Finger,  M. ,  and  Kurrle,  J.  E.,  "Equation  of 
State  of  Detonation  Products,"  Proc.  5th  Symposium  (Int.)  on  Detonation,  ONR 
ACR-184  (Washington,  D.  C.:  U.S.  Gov.  Print.  Office,  19/2),  pp.  503-5l2. 


FIGURE  1  DETONATION  BEHAVIOR  OF  HOT-FRES5ED  TNT  (REF  8) 
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which  flow  easily  under  pressure.  ftempel6  has  demonstrated  the  effect  for 
TNT  pressed  at  72-76°C  In  the  presence  of  a  llttl®  acetone.  Figure  1 
Illustrates  It  for  8-13  mm  dla  charges.  Since  Fig.  1  shows  the  trend  dc 
Increasing  with  increasing  p0,  the  opposite  of  that  found  for  cold  pressed 
charges'.  It  Is  reasonable  to  suppose  that  the  hot-pressing  (aided  by  solvent) 
has  not  only  Increased  but  also  considerably  changed  the  physical  nature 
of  the  charge.  This  Is  particularly  likely  In  view  of  the  low  melting  point 
of  TNT  (80.5°C).  In  other  words  the  very  low  porosity.,  hot-pressed  charges 
approach  a  cast  TNT  rather  than  a  cold  pressed  one  In  their  structure.  While 
It  has  not  been  established  that  cast  and  pressed  TNT  differ  In  0,  It  has 
been  reported  that  at  *30  mm  dla  and  p0  *  1.62  g/cm3  the  cast  TNT  has  a 


Gibbs,  T.  R. ,  and  Poolato.  A..  eds..  LASL  ExdIosIvp  Prnnertv 


(Berkeley:  U.  Cal.  Press,  1980). 
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SINGLE  EXPLOSIVES 


TNT 

Ref.  4  reports  what  are  probably  the  most  precise  detonation  volocltles 
measured  for  1NT.  However,  all  charges  of  p«  >  1.44  g/cnr  were  made  by 
hot-pressing  (,'70oC).  The  equation  derived  for  Infinite  diameter  or  Ideal 
detonation  velocity  In  the  range  p0  *  0.9  g/cm3  was 


1.8727  +  3.1872  pQ  0.9  s  pQ  s  1.5342 


6.7625  +  3.1872  (pQ  -  1.5342)  (1) 

-  25.102  (p0  -.  1.5342)2 

+  115.056  (p0  -  1.5342) 3  1.534*  s  p0  i  1.636  g/cm3 

q  in 

Recent  handbooks  *  quote  Eq.  1  without  the  cubic  term  or  any  comment  on  Its 
absence. 

Fig.  2  compares  Eq.  1  with  that  derived  by  Hurwitz:^ 

0.  »  1.785  +  3.225  pQ  (2) 

Each  curve  has  been  terminated  at  the  highest  experimental  density  reported. 

The  two  sets  of  results  are  within  experimental  error  except  at  the  higher 
densities  where  hot-pressed  charges  were  used.  There  Eq.  1  indicates  a 
sharp  decrease  In  the  slope;  that  could  result  from  an  increased  critical 
diameter  such  as  shown  in  RempseVs  data.  An  increased  would  Increase 
the  diameter  effect  on  the  measured  detonation  velocity  at  any  given  diameter 
and  hence  might  result  In  D(neasured  <  Dr 

10Qobratz,  B.  M. ,  ed.,  LLNL  Explosives  Handbook.  UCRL-S2997,  1981  with  1982 
addenda. 

^Campbell,  A.  W. ,  and  Engelke,  R. ,  "The  Diameter  Effect  in  High  Density 
Heterogeneous  Explosives,"  Sixtn  Symposium  (Int.)  on  Detonation,  ONR  ACR  221, 
(Washington,  D.  C.:  U.  S.  Gov.  Print.  Office,  1978)  pp.  642-663. 


D^ (mm/us)  « 
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FIGURE  2  0(^Q)  CURVES  AND  DATA  FOR  TNT 
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Ref.  9,  In  addition  to  Eq.  1  without  the  cubic  term,  also  quotes 
7.045  mm/us* as  Dj  for  Dressed  TNT  as  p0  ■  1.620  g/cm3.  However  Eq.  1 
produces  0^ (1.620  fl/cm3)  •  6.924  mm/us  (6.051  mm/us  if  cubic  tern  omitted). 

If  the  change  in  slope  shown  by  Eq.  1  is  Ignored,  the  relation  becomes 

Di  ■  1.8727  ♦  3.1872  p0  pQ  i  0.9  g/cm3  (3) 

Moreover,  for  p0  ■  1.62  g/cm3,  Eq.  3  produces  D^  •  7.036  mn/ws,  In  quite  good 
agreement  with  the  7.045  rnn/us  value.  Typical  high  density  experimental  values 
of  (Oi,p0)  are:  (7.00,1.62)8,  (6.91,1.59)“,  and  (6.91,1.64) ;13  the  first  two 
sets  are  Russian,  the  last  from  LANL  (hot-pressed  charge).  These  points  are  also 
plotted  In  Fig.  2  to  demonstrate  the  pattern  of  Russian  data  falling  to  exhibit 
the  change  of  slope  at  high  density.  NSWC  has  obtained  TNT  at  1.60  g/cm3  by 
pressing  cold  charges  In  the  hydraulic  press  and  at  1.64  g/cm3  in  the  Isostatic 
press.  It  seems  possible  and  even  likely  that  the  Russian  data8*12  were 
obtained  by  pressing  cold  charges.  In  view  of  the  possible  effect  of  hot- 
pressing  on  the  charge  performance.  It  Is  suggested  that  Eq.  3 
be  used  for  cold-pressed  charges.  It  will  be  so  used  with  the  prediction 
method  for  mixtures  described  later.  Of  course,  an  experimental  Investigation 
of  the  measurable  effects  on  dr  and  D(d)  of  hot-pressed  vs  cold-pressed  charges 
is  strongly  recommended  as  well. 


PETN 


PETN  has  a  higher  melting  point  than  TNT  (141°C),  but  It  too  flows  easily 
under  pressure.  The  most  precise  0  measurements  on  PETN  are  probably  those 
of  Ref.  5.  That  reference  did  not  state  whether  the  pressing  was  of  cold 
or  heated  PETN.  However,  LLNL  generally  uses  hot-pressing  as  does  LANL. 


Ref.  5  presents  the  vs  p0  curve  as 


(mm/us)  « 
* 


2.14  +  2.84  p0 
3.19  +  3.70  (p( 
7.92  +  3.05  (p, 


0.37) 

1.65) 


p0  <  0.37  g/cm* 

0.37  <  p0  <  1.65  (4) 

P0  >  1.65 


12 

Dremln,  A.  N.,  and  Shvedov,  K.  K. ,  "Determination  of  C-J  Pressure  and 
Reaction  Duration  In  the  Detonation  Wave  of  HE,"  PMTF  No.  2,  1964, 
pp.  154-159,  through  a  translation. 

13Mader,  C.  L.,  Los  Alamos  Rept.  LA-2900,  1963. 

♦This  value  was  derived  from  Ref.  8  after  smoothing  the  function  C(d)  according 
to  the  method  of  Ref.  11. 
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In  the  Intermediate  density  range  of  1.0  t  0.2  g/cm  ,  Eq.  4  and  Refs.  1, 

14.  and  IS  are  In  very  close  agreement.  Ref.  1  limits  Its  data  to  o0  a  0.8 
g/cm.  Below  that  limit.  Refs.  14  and  IS  are  still  In  good  agreement  with  .. 
Eq.  4  provided  their  data  are  forced  to  extrapolate  to  the  theoretical  value10 
of  2.19  mm/»s  at  p0  •  0.01  g/cm3;  that  value  compares  well  with  2.17  mm/us 
at  p0  ■  0  from  Eq.  4.  It  was  probably  the  difficulty  of  obtaining  uniform 
and  stable  charges  at  p0  <  O.S  g/cc  that  resulted  In  reported  D  measurements 
at  very  low  p0  well  above  the  Eq.  4  curve  In  Ref.  14. 

3 

In  the  range  p0  a  0.95  g/cm  ,  Eq.  4  Is  shown  as  the  dashed  line  of 
Fig.  3.  Also  shown  Is  the  curve 

01  •  1.600  +  3.950  pQ  (5) 

from  Ref.  1.  Ref.  5  quotes  six  experimental  values  attributed  to  a  report 
by  E.  A.  Christian  and  H.  G.  Snay.  However,  those  values  were  all  derived 
from  Eq.  5.  Table  1  contains  the  additional  experimental  data  plotted 
In  Fig.  3.  They  are  In  good  agreement  with  Eq.  4  up  to  p0  »  1.65  g/cm3. 

Above  that  density,  LANL  data  lie  below  and  the  Russian  data  lie  above  the 
Eq.  4  curve.*  Foi  the  same  reasons  considered  In  the  case  of  TNT,  It  Is 
suggested  that  the  equation 

D1  «  1.82  ♦  3.70  p0  pQ  >  0.8  g/cm3  (6) 

be  used  for  computing  0  of  mixtures  containing  PETN. 

Other  Pure  HE 


Other  single  explosives  to  be  considered  here  are  8TNEU,  DINA,  RDX(HMX), 
tetryl  and  TNETB. 

BTNEU,  bls(trinitroetnyl)  urea,  O-CKNHCHgCfNOgh^  Is  a  high  energy, 
sensitive  explosive  with  a  voldless  density  of  1.86  gran3  and  a  decomposition 

14Fr1eder1ch,  W.,  ,;Uber  die  Detonation  der  Sprengstoffe,"  Z.  ges.  Schless- 
Sprengstoffw. .  Vol .  28,  1933,  pp.  2-6,  51-3,  80-3  and  113-18. 

15 

Paterson,  S. ,  Teknlck-Vetenskaplig  Forskning,  Vol.  29,  1958,  pp  109-120. 

The  same  data  are  quoted  In  the  text,  0.  Taylor,  Detonation  in  Condensed 
Explosives  (Oxford:  Clarendon  Press,  19501. 

^Stesik,  L.  N.,  and  Shvedova.,  N.  S.,  "Detonation  of  Condensed  Explosives 
with  Low  Charge  Densities,"  PMTF,  No.  4,  1964,  pp.  124-126,  ONI  translation 
No.  2111. 

*It  is  of  interest  that  the  LANL  curve  for  PETN® 

D  »  1.608  +  3.933  pQ  0.57  <  pQ  <  1.585 

extrapolates  to  8.4Q  nrn/ys  at  p  *  1.75  g/cm3.  This  value  lies  quite  close 
to  the  Russian  data.  Of  course,  this  is  also  true  of  Eq.  6- 


FIGURE  3  D  (^0)  CURVES  AND  DATA  FOR  PETN 
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TABLE  1 

Experimental  Values  for  D^(pQ)  of  PETN  in  Range  pQ  >  0.95  g/cm 


‘Jo  3 

g/cm 

u 

mm/us 

Ref. 

1.77 

8.60 

17 

1.73 

8.35 

14 

1.67 

7.98 

13 

1.66 

8.10 

12 

1.65 

7.92 

18 

1.62 

7.91 

14 

1.51 

7.42 

12 

1.51 

7.42 

14 

1.37  . 

6.97 

14 

1.03 

5.62 

•  14 

0.97 

5.33 

14 

0.95 

5.30 

12 

*7Apin,  A.  Ya. ,  and  Voskoboynikov. ,  I.  M. ,  "Calculation  of  Parameters  of  a 
Detonation  Wave  of  Condensed  Explosives,"  PMTF,  No.  4,1960,  pp.  54-55,  ONI 
Translation  No.  2083.  * 

18Craig,  B.,  through  Ref.  5. 
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(7) 


point  of  o-186°C.  Its  D-j  was  derived  from  experimental  measurements  in 
Refs.  19  and  20: 

Di  =  1.42  +  4.'08p 

The  charges  used  were  highly  porous  or  contained  2%  wax  as  a  binder  at 
^1.6  g/cm3.  However  Eq.  7  extrapolates  to  within  0.1%  of  the  value  measured 
at  1.86  g/cm3  in  small  scale  (confined)  charges. 

DIMA,  diethanolnitramine  dinitrate,  OgN-N ( CH2CH20N02 ) 2  bss  a  voidless 
density  of  1.67  g/cm3  and  m.p.  of  52.5°C.  It  is  an  easily  synthesized, 
easily  castable  energetic  explosive,  and  has  been  used  in  DDT  studies. 

Ref.  1  gives  for  this  HE 

Di  *  3.00  +  2.95p0  (8) 

A  review  has  been  made  of  the  data  available  since  Ref.  1  was  published. 

These  data  are  collected  in  Table  2  and  plotted  in  Fig.  4.  Measured  D  for 
cast  DINA  was  Included  because  Ref.  1  showed  that  its  measured  velocity 
fell  on  the  curve  obtained  with  pressed  charges. 

The  solid  line  of  Fig.  4  is  a  least  squares  fit  to  the  eight  data 
pairs  of  Table  2.  Its  equation  Is 

Di  *  3.03  +  2.93p0  (9) 

with  a  *  0.66  mm/us.  There  is  consequently  no  justification  for  changing 
Eq.  8,  particularly  since  it  treats  a  greater  amount  of  data  from  charges 
prepared  by  standardized  methods  at  ERL. 

RDX(HMX)  share  the  same  D(pQ)  function  which  is 

D.  *  2,56  +  3.47p0  (10) 

from  Ref.  10.  This  equation  fits  the  RDX  mixtures  data  at  p0  »  1.0  g/cm3 
better  than  that  in  Ref.  9  although  both  give  equivalent  values  at  voidless 
density.  It  was  demonstrated  some  time  ago  that  the  low  porosity  LANL  data9 


19 

August,  W.  T.,  Liddiard,  T.  P.,  Lovenberg,  C.  C.,  Shapiro,  N.,  Solem,  A.  D., 
Stresau,  R.  H.,  and  Svadeba,  G.,  "Evaluation  of  BTNEU  as  a  Substitute  for  RDX," 
NavOrd  2287,  Dec  1951. 

20 

Coleburn,  N.  L.,  and  Liddiard,  T.  P.,  "The  Rates  of  Detonation  of  Several 
Pure  and  Mixed  Explosives,"  NavOrd  2611,  Sep  1952. 
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TABLE  2 

(pQ)  Data  Recently  Reported  for  DINA 


/  D 
mm/us 

Ref. 

1.67 

8.00 

17 

1.64 

7.80 

21 

1.60  c* 

7.72 

22 

1.60 

7.73 

23 

1.55 

7.58 

23 

1.48 

7.40 

17 

1.36 

7.00 

21 

0.95 

5.80** 

21 

9 1 

Dremln,  A.  N.,  Savrov,  S.  D. ,  Trofimov,  V.  S. ,  and  Shvedov,  K.  K.,  Detonation 
Maves^  in  ^Condensed  Media  (Moscow:  Izd-vo  Nauka,  1970)  FTD-HT-23-16M-71 

22 

Price,  D. ,  Jaffe,  I.,  and  Toscano,  J.  P. ,  “Development  of  the  Continuous 
Wire  Method  II,"N0LTR  66-21,  Mar  1966. 

23 

Fedoroff,  B.  T.»  and  Sheffield,  0.  E. ,  eds.,  Encyclopedia  of  Explosives  and 
Related  Items.  Vol.  5  (Dover,  N.  J.:  PI catinny  Arsenal ,1^2),  p.  D-l24i. 

♦Only  charge  for  which  preparation  data  are  available. 

♦♦Value  given  in  translation  used  Is  5.08.  However,  this  is  clearly  a 
typographical  error  as  can  be  checked  with  values  of  detonation  pressure 
and  particle  velocity  that  are  also  listed. 
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(and  hence  the  low  porosity  ILNL  data^)  for  ROX  were  more  consistent  than 
those  of  NSWC?* 

Tetryl  Is  another  Important  commercial  explosive  for  which  the  0(po) 
function  does  not  seem  to  have  been  reexamined  since  19451.  The  relation 
given  In  Ref.  1  Is 


D1  «  2.375  ♦  3.225pq  (11) 

On  the  other  hand,  Ref.  9  gives 

0^  »  2.742  +.2.935p0  1.3  s  P0  S  1.69  (12) 

which  Is  appreciably  different  from  Eq.  11.  Table  3  contains  some  of  the 
more  recent  data  for  tetryl,  and  Fig.  5  displays  them  together  with  Eqs.  11 
and  12.  It  Is  evident  that  the  data  conform  better  to  Eq.  12  than  to  Eq.  11; 
hence  the  former  will  be  used  for  mixtures.  Ref.  21  values  seem  to  run  a 
bit  low  and  Ref.  17,  a  bit  high  as  they  did  for  DINA.  Aside  from  that,  the 
larger  discrepancies  at  low  porosity  might  arise  from  different  methods  of 
charge  preparation,  different  Instrumentation,  or  both. 

TNETB,  trlnltroethyltrlnltrobutyrate,  (NOjtaCC^C^COOCHj^NOzh,  Is 
a  castable  (m.p.  93°C)  high  energy  explosive.  Since  Its  D(pp)  function, 
as  reported  In  the  1965  NSWC  Explosives  Handbook,  Is  somewhat  In  error, 
the  corrected  function  Is  presented  here  for  convenience;  Its  derivation 
was  reported  In  TR  68-182. 

Di  »  1.947  +  3.660po,  <r  -  0.015  rnn/us  (  (13) 

Eq.  13  Is  very  similar  to  Eq.  6,  the  analogous  relation  for  PETN. 


24 

Price,  0.,  "Inter-relationship  of  Explosive  Characteristics  III"  NavOrd 
4510,  Apr  1957. 
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TABLE  3 

Addition*!  D^(oQ)  Data  for  Tetryl 


mm)us 

Ref. 

1.71 

7.85 

10 

1.70 

*  7.860 

17 

1.70 

7.560 

13 

1.68 

7.50 

21 

1.614 

7.581 

25 

1.60 

7.400 

23 

1.55 

7.30 

23 

1.51 

7.17 

26 

1.506 

7.150 

27 

1.44 

6.875 

28 

1.36 

6.68 

21 

1.22 

6.291 

14 

0.9-0.95 

5.36 

21 

0.95 

5.390 

28 

05 

Colebum,  N.  L. ,  "C-0  Pressures  of  Several  Pure  and  Mixed  Explosives," 
NOLTR  64-58,  Sep  1964. 

26 

Clalrmont,  A.  R. ,  and  Oaffe,  I.,  "Analysis  of  the  Optical  Determination  of 
Detonation  Velocity  In  Short  Charges,"  NOLTR  64-23,  May  1964. 

27 

Cybulskl,  W.  B. ,  Payman,  W. ,  and  Woodhead,  D.  W. ,  Roy.  Soc.  London  Proc.. 
Vol.  197A,  1949,  pp.  51-72.  — 

£®Woodhead,  D.  W. ,  and  Tltman,  H. ,  "Detonation  Phenomena  In  a  Tubular  Charge 
of  Explosive,"  Exploslvstoffe,  Vol.  5,  113-123  and  6,  141-155,  1965. 

29Pr1ce,  D.,  Clalrmont,  Jr.,  A.  R.,  Erkman,  J.  0.,  and  Edwards,  D.  J., 

"Ideal  Detonation  Velocity  of  Ammonium  Perchlorate  and  Its  Mixtures  with 
HE,"  NOLTR  68-182,  Dec  1968. 
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FIGURE  8  D  ( pQ)  CURVES  ANO  DATA  FOR  TETRYL 


a  y  a  a 


MIXTURES  OF  HE 


The  experimental  values  of  the  mixed  HE  will  be  obtained  from  the  curves 
reported  by  Colebum  and  Llddlard20.  In  general,  deriving  such  curves  results 
In  averaging  and  compensating  for  errors  Introduced  by  unrecognized  differences 
In  charge  preparation,  measurements,  and  record  reading.  Thus  a  velocity 
value  on  a  D-pg  curve  derived  from  a  number  of  experimental  values  over  the 
density  range  Is  considered  of  equal  or  greater  accuracy  than  the  Individual 
experimental  values  upon  which  It  Is  based.  (Similarly  en  average  of  two 
or  more  experimental  measurements  on  charges  of  the  same  density  Is  considered 
to  be  more  accurate  than  Individual  measurements.)  This  Is  also  considered  to  be 
.the  case  for  pure  explosives. 

Table  4  contains  a  summary  of  the  D(p0)  functions  chosen  for  the  single 
HE  and  those  given  In  Ref.  24  for  the  mixtures.  All  are  In  the  form 

Di  *  a  ♦  bpQ  '  (14) 

9 

In  addition  to  a  and  b,  the  Intercept  and  the  slope,  pc  (the  voldless  or 
theoretical  maximum  density  or  TMO)  Is  also  listed.  Eq.  14  can  be 
rearranged  to 

Di  ■  a  +  bp_A  (15) 

where 

A  -  pQ/pc  -  XTMD/100 

For  this  reason,  bpc  is  also  listed  In  Table  4. 

In  addition.  Table  4  contains  the  curves  derived  from  the  data  of 
Ref.  20  for  mixtures  of  HE.  Only  three  of  the  mixtures  are  Included  In 
the  Ref.  1  compendium:  pentollte,  cyclotol  60/40  (Ref.  1  does  not 
distinguish  between  this  cyclotol  and  Comp  B)  and  tetrytol  70.  In  the 
first  two  cases,  the  curves  were  essentially  the  same  l.e.,  less  than  IX 
difference  for  pentollte  and  l  -  1.4X  difference  for  the  cyclotol;  the 
Ref.  20  data  gave  higher  values,  possibly  because  Comp  B  (Ref.  1)  Includes 
IX  wax.  However,  the  two  tetrytol  curves  differed  significantly  (ca.  5X 
at  1  g/cm3  -  not  as  much  at  higher  densities.)  The  curve  from  Ref.  20 
was  used. 
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To  predict  0^  for  a  mixture  from  known  D^(p0)  of  the  components,  the 
additivity  relation 

V’oW  “  J  *j  .  «W8 

where  xj  Is  weight  fraction  of  component  j  and  Di  Is  the  Ideal  detonation 
velocity  of  consonant  j  at  the  porosity  of  the  mixture,  was  suggested  some 
time  ago28.  It  was  at  that  time  applied  only  to  materials  at  voldless 
density;  for  the  present  application  It  has  been  generalized  to  any  porosity 
In  the  linear  ranee  of  the  D(p0)  functions.  Eq.  1$  In  conjunction  with  Eq.  IS 
produces  the  predictive  equation 

°1(po>mix  *  2  *j  |\  *  bj  Pcjp)  ^  J  (17) 

Eq.  17  has  been  used  to  predict  the  detonation  velocity  of  mixtures  at 
densities  of  1.0  and  1.6  g/ern3.  The  predicted  and  observed  values  are 
listed  In  fable  4. 

About  forty  years  ago,  Urlzar  formulated  an  empirical  equation  for 
predicting  detonation  velocity  from  a  knowledge  of  detonation  velocity  at 
crystal  density  (Dc)  of  each  mixture  component.  His  expression,  described 
on  p.  8*10  of  Ref.  10  Is 

W  •  J  yJ  Dcj  (l8) 

where  yj  Is  volume  fraction  of  component  j.  For  non*explos1ve  components 
such  as  air,  A1  and  NaCl,  he  gave  empirical  values  for  Dc.  For  air  (voids) 
this  empirical  number  Is  1.5  mm/ws. 

Since  for  each  solid  component. 


and  for  the  voids 


solids 

y#  •  1  -  I  yj  »  1  *  (*  TMD/100) 


(19) 

(20) 


Eq.  18  becomes 

w  ■  •«  Ts  %  (*j  *  bi  pcj> +  »•* ».  <n> 

Eq.  21  was  used  to  calculate  Dj  at  p0  ■  1.0  and  1.6  g/cm3  for  both  single 
explosives  and  mixtures.  The  experimental  and  computed  values  are  listed 
In  Table  5. 
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The  last  approximation,  again  applicable  to  either  pure  HE  or  their 
mixtures,  to  be  considered  here  Is  that  of  Kamlet30.  This  approximation 
1$  based  on  the  arbitrary  decomposition  that  uses  the  available  oxygen 
In  the  HE  to  form  H^O(g)  and  COg  In  that  sequence.  Then 

D1  (mm/us)  «  1.01*1/2  (1  +  1.30  p0)  (22) 

where  p0  Is  density  In  g/cm3,  and  *  ■  NM1^2  qJ£?;  N,  M,  and  Q*rb  are, 
respectively,  number  of  moles  of  gaseous  prodflcts  per  gram  HE,  average 
molecular  weight  of  gaseous  products,  and  the  specific  chemical  energy  of 
the  detonation  reaction.  The  values  of  N,  M,  Qarb»  and  *  are  all  determined 
by  the  arbitrary  equation  used.  Again  values  at  p0  ■  1.0  and  1.6  g/cm3 
have  been  computed  and  are  compared  to  the  observed  values  In  Table  6. 

The  Table  4  data  show  that  the  additivity  method  of  Eq.  17  Is  a  very 
good  way  to  predict  velocities  of  these  mixed  HE  at  moderate  to  low  porosities 
(pQ  *  1.6  g/cm3).  Its  average  absolute  error  Is  less  than  IX  and  the 
Individual  error  less  than  2%.  With  the  exception  of  mixtures  containing 
PETN,  it  Is  equally  effective  at  high  porosity  (p  *  1.0  g/cm3),  but  the 
PETN  mixtures  show  .errors  of  up  to  -3.5%.  A  possible  reason  for  this  Is 
that  Eq.  17  assumes  no  Interaction  between  the  detonation  products  of  the 
low  and  high  oxygen  components  of  the  mixtures.  At  high  porosities  the 
reaction  time  Is  greater  and  the  reaction  temperature  higher  than  at  low 
porosity.  Under  these  conditions  Interactions,  absent  at  low  porosities, 
might  occur  and  result  In  Increasing  the  detonation  velocity.  This 
possibility  seems  somewhat  probable  In  view  of  the  fact  that  the  Kamlet  method, 
which  assumes  complete  Interaction,  Is  satisfactory  for  these  mixtures  at 
Po  *  1*0  g/cm3  despite  Its  unsatisfactory  result  at  higher  density  for  the 
high  oxygen  component,  PETN30. 

As  the  Table  5  data  show,  the  Urlzar  method  is  almost  equally  acceptable 
for  both  pure  and  mixed  HE  at  pp  %  1.6  g/cm3.  It  Is,  however,  quite 
Inadequate  at  high  porosities  where  Individual  results  vary  from  the  experi¬ 
mental  by  as  much  as  10%  and  the  average  absolute  difference  Is  6-7%. 

The  Kamlet  method  (Table  6)  gave  excellent  results  for  this  set  of 
mixtures  at  p0  ■  1.6  g/cm,  and,  with  the  exception  of  the  60/40  BTNEU/TNT, 
acceptable  ones  at  the  higher  porosity.  However,  Its  prediction  for 
pure  explosives  was  not  good  at  the  lower  density.  With  the  exception 
of  DINA,  the  HE  showing  the  higher  discrepancies  at  p0  ■  1.6  g/cm3  exhibited 
greater  discrepancies  (up  to  10%)  at  the  greater  porosity.*  Those  HE  showing 


30Kamlet,  M.  0.,  and  Hurwltz,  H.,  J.  Chem.  Phys.,  Vol.  No.  48,  No.  8,  1968, 
pp.  3685-3692. 

♦Comparisons  In  Tablgs  4  and  5  are  at  fixed  densities.  If  Instead  fixed 
porosity  is  used,  the  picture  is  the  same.  E.g.,  at  57%  TMD,  instead  of 
p0  ■  1  g/cm3,  the  single  HE  show  an  average  absolute  error  of  5.5  and 
3.6%,  respectively  for  values  from  the  Urlzar  and  Kamlet  methods. 
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the  greatest  differences  between  predicted  and  observed  D  were  BTNEU,  PEfN, 
and  TNETB  -  all  high  oxygen  HE.  The  discrepancy  for  PETN  was  noted  in 
Ref.  30,  but  no  data  for  high  porosity  BTNEU  and  TNETB  were  considered.* 

Of  course,  results  of  the  comparison  between  predictions  by  this  method 
and  the  experimental  values  will  be  influenced  by  the  choice  made  for  the 
experimental  curve  more  than  the  two  earlier  comparisons  because  the  other 
methods  are  both  based  on  experimental  values  whereas  this  one  is  not. 


*In  the  case  of  DINA,  the  experimental  data  used  in  Ref.  30  were  the  Ref.  17 
data  shown  on  Fig.  4.  The  higher  density  value  was  well  above  the  curve 
chosen  here  for  the  linear  curve;  hence  the  discrepancy  observed  in  Ref.  30  was 
not  as  great  as  the  present  +3.1%. 
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SUMMARY 


It  is  suggested  that  the  recently  reported  sharp  changes  in  slope  of  the 
D(p0)  curves  at  nigher  densities  for  TNT*  and  PETN5  may  result  from  changes 
in  the  charges  caused  by  hot-pressing.  In  other  words,  the  increased 
homogeneity  of  the  hot-pressed  charges  can  Increase  the  critical  diameter 
and  thereby  the  diameter  effect  (a  lowering  of  D  at  finite  diameters)  on 
measured  detonation  velocity.  Therefore,  it  is  strongly  recommended  that 
the  effect  of  hot-pressing  on  0  be  investigated. 

D(p0)  data  for  TNT,  PETN,  BTNEU,  DINA,  RDX(HMX)  and  tetryl  have  been 
reviewed  and  a  best  experimental  curve  selected  In  each  case.  These  curves 
have  been  used  to  predict  the  D  values  of  HE  mixtures  by  an  additivity  rule. 
The  predicted  values  were  In  excellent  agreement  with  the  measured  values 
(Ref.  20)  at  high  and  low  porosities  except  for  two  high  porosity  mixtures 
containing  PETN.  The  underestimate  In  these  two  cases  might  result  from 
interaction  of  the  detonation  products  of  PETN  with  those  of  lower  oxygen 
components. 

The  same  single  HE  data  were  used  to  obtain  Dc  values.  These  in  turn 
were  used  In  Urlzar's  approximation  for  D-j.  It  gave  very  good  results  at 
p0  *  1*6  g/cm3;  very  poor  ones  at  p0  *  1.0  g/cm3. 

Kamlet's  approximation  for  D-j  from  the  H2O  -  CO?  arbitrary  decomposition 
produced  values  for  the  mixtures  in  good  agreement  with  experiment  at 
p0  *  1.6  g/cm3;  fair  agreement  at  p0  *  1.0  g/cm3.  However,  the  same, 
procedure  resulted  in  over  three  percent  difference  at  pg  *  1.6  g/cnr  for 
two  single  explosives  and  in  large  differences  (up  to  10%)  for  four  single 
explosives  at  p0  «  1.0  g/cm3.  The  largest  differences  were  for  three  high 
oxygen  content  HE:  BTNEU,  PETN,  and  TNETB. 

The  best  overall  method  for  predicting  D  of  mixtures,  satisfactory  except 
for  high  porosity  PETN  mixtures.  Is  the  additivity  method  of  Eq.  17.  Its 
failure  in  this  case  is  attributed  to  reaction  between  the  detonation  products 
of  PETN  and  those  of  oxygen  deficient  components. 
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